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ABSTRACT 

The eccentricity distribution of exoplanets is known from radial velocity surveys 
to be divergent from circular orbits beyond 0.1 AU. This is particularly the case for 
large planets where the radial velocity technique is most sensitive. The eccentricity of 
planetary orbits can have a large effect on the transit probability and subsequently 
the planet yield of transit surveys. The Kepler mission is the first transit survey that 
probes deep enough into period-space to allow this effect to be seen via the variation 
in transit durations. We use the Kepler planet candidates to show that the eccentricity 
distribution is consistent with that found from radial velocity surveys to a high degree 
of confidence. We further show that the mean eccentricity of the Kepler candidates 
decreases with decreasing planet size indicating that smaller planets are preferentially 
found in low-eccentricity orbits. 
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1 INTRODUCTION 

Planets discovered using the radial velocity (RV) method 
have dominated the total exoplanet count until recently, 
when the transit method has made increasing contributions. 
The long time baseline of RV surveys has allowed the de- 
tection more diverse orbital geometries than achievable by 
ground-based transit surveys. The Kepler mission, however, 
with its multi-year baseline, can begin to probe into parame- 
ter space previously reserved for RV studies. At longer peri- 
ods, orbits tend to diverge significantly fr om the circular cas e 
beyond a semi-major axis of ~ 0.1 AU (jButler et al.ll2006r ). 
although there m ay be small observational biases that skew 
this distribution (|Shen fc Turner 2008 ). This insight has led 
to numerous attempts to account for eccentric ity in the con- 
text of planet formation and orbital stabil ity (Ford fc Rasid 



2008; Malmberg & Davics 2008; Matsumura et all 200. 



Juric fc Tremainel 120081 : IWang fc FordllffllH') and the influ 



ence of tidal circularization ( Pont et al.ll2011^ 



It has been shown h ow eccentricity distribut i on ef - 
fects transit probabilities (jKane fc von Br aun 2008, 2009) 
and projected yields of transit surveys (|Barnea.2007. : .Burka 



120081 ). This influence is minor for the ground-based sur- 
veys since they are primarily sensitive to giant planets in 
short-period orbits. However, the Kepler mission is expected 
to be impacted by this distribution since it probes out to 
much longer periods with a much reduced disadvantage 
of a wi ndow function that aff ects observations from the 
ground (Ivon Braun et all |2009| ) . A comparison of the Ke- 



pler results in the context of eccentricity and transit du- 
rations with the RV distribution has been suggested by 
I Ford e t aU (^2008) and 'Zaka rnska et al.l (|201ll ) and carried 
out by Moorhead ct al. (2011J), but initial planet candidate 
releases by the Kepler project do no t provide enough pe- 
riod sensitivity (jBorucki et al.l l2011alrb|). The most recent 
release of Kepler planet candidates bv lBatalha et all (|2012l ) 
increases the total number of candidates to more than 2,300 
and the time baseline probed to beyond 560 days. This 
has several implications for studies of eccentricity distribu- 
tions. The Kepler mission is sensitive to planets significantly 
smaller than those accessible by current RV experiments and 
thus allows a more in-depth study of the dependence of ec- 
centricity on the planet mass/size and multiplicity. If the ec- 
centricity distributions of Kepler and RV planets were found 
to be substantially different then this may reveal a selection 
effect in the way Kepler candidates are selected which is 
biased against eccentric orbits. A direct comparison of the 
two distributions, provided they are consistent for the planet 
mass/size region where their sensitivities overlap, will allow 
a more exhaustive investigation of orbital eccentricity to be 
undertaken. 

Here we present a study of the eccentricity distribution 
of planets discovered with the RV method and the complete 
list of Kepler planet candidates. We calculate expected tran- 
sit durations for circular orbits and compare them with ei- 
ther calculated or measured eccentric transit durations (32)). 
Our results show that the measured transit durations from 
RV data (!j3]) and the Kepler candidates (!j4]) are consistent 
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with having the same distribution. We estimate the impact 
parameter distribution for the Kepler candidates and show 
that their mean eccentricity decreases with decreasing planet 
size ([J5|, which supports the hypothesis that smaller planets 
tend to be found in multiple systems in near-circular orbits. 
We discuss additional astrophysical aspects in [JS] and con- 
clude in ijTl 



2 ECCENTRICITY AND TRANSIT 
DURATION 

A concise description of exoplanetary transit model- 
ing and associated parameters is presented else where 
l|Mandel fc Agolll2002l : ISeager fc Mallen-Ornelas "2002). Here 
we concentrate on the relevant details to our analysis: transit 
duration and eccentricity. The transit duration for a circular 
orbit is given by 



P / \/('R* + -Rp)^-a^cos2i 

tcirc ~ — arcsin — 

TT I a 



(1) 



where P is the orbital period, a is the semi-major axis, i is 
the orbital inclination, and _R* and Rp are the stellar and 
planetary radii respectively. The impact parameter of a tran- 
sit is given by 

b = — — cos i (2) 

and is defined as the projected separation of the planet and 
star centers at the point of mid-transit. 

For an eccentric orbit, the star-planet separation r is 
time-dependent and is given by 

1-fecos/ ^^' 

where e is the orbital eccentricity and / is the true anomaly. 
Replacing a with r at the time of inferior conjunction in 
equations [l] and [2] provides generalized expressions for tran- 
sit d uration and impact parameter for non-circular orbits. 
iBurkG (,2008,) converts these expressions into the scaling fac- 
tor 

tecc _ Vl — e^ 

t^ ^ l + ecos{u}- 90°)' 

where w is the periastron argument of the orbit. 



(4) 



3 RADIAL VELOCITY ECCENTRICITY 
DISTRIBUTION 

We first investigate the eccentricity distribution of the plan- 
ets discovered with the RV technique and the subsequent im- 
pact on the predicted transit duration. The Exoplanet Data 
Explorer (EDEJj stores information only for those planets 
that have co mplete orbital soluti ons and thus are well suited 
to this study (| Wright et aLlbOllI ). The EDE data are current 
as of 2012 February 24 and include 204 planets after the fol- 
lowing criteria are applied: logg > 3.5 to exclude giant stars 
and a < 1.5 AU to produce a sample that covers the same 
region in parameter space as the Kepler candidates. 



To calculate the transit duration, one needs an estimate 
of the planetary radius. For planets that are not known to 
transit, we approximate the planetary r a dius u sing the sim- 
ple model described bv iKane fc Gelind (120121 ). This model 
adopts a radius of 1 Jupiter radius for masses > O.SMjupUer 
and utilizes a power law fit to the masses and radii of the 
known transiting planets for masses < O.SMjupiter- In order 
to estimate the radius of the host star, we use the following 
relation related to the surface gravity 



logs = log 



+ log 5© 



(5) 



where logg© = 4.4374 (jSmallevll2005l ). Using the equations 
of Section [2] and the measured orbital parameters, we calcu- 
late the transit duration for both the circular and eccentric 
cases. We then take the absolute value of the difference be- 
tween the two durations as a diagnostic for the eccentricity 
distribution. 

The top-left panel of Figure[T]shows the eccentricity dis- 
tribution of the RV planets taken from EDE as a function 
of a. The distribution begins to diverge from mostly circular 
orbits beyond 0.04 AU, and by 0.1 AU it has an eccentric- 
ity range of 0.0-0.5. This distribution is widely attributed 
to tidal damping of the orbits after the disk has dissipated. 
According to Goldreich fc Sotcr (1966), the timescale for or- 
bital circularization is oc a^'^M^^'^ where M* is the stel- 
lar mass. Note that the two planets inside of 0.04 AU with 
e > 0.1 are GJ 436b and GJ 581e, both of which have M 
dwarf host stars. 

The other three panels in Figure [T] show the calculated 
duration difference by At as a function of a. At is the ab- 
solute value of the difference between the calculated tran- 
sit duration for a circular orbit and the calculated transit 
duration based upon the measured orbital parameters (i.e.. 
At = |tcirc — iecci), which is indicative of the divergence from 
the assumption of only circular orbits. The top-right panel 
assumes edge-on orbits (i = 90° ; 6 = 0) for both the circu- 
lar and eccentric cases. We show the effect of increasing the 
impact parameter of the transits for & = 0.5 and b = 0.8. 
The mean of At is not significantly changed except for rel- 
atively high values of b. We evaluate the significance of this 
distribution in the following section. 

It should be noted that, in order for a transit to take 
place, i can at most only slightly deviate from 90°. The con- 
sequent small angle approximation means that the uniform 
distribution of i values maps to a uniform distribution of b 
values, making all values of b equally likely to occur. 



4 ANALYSIS OF KEPLER TRANSIT 
DURATIONS 

The release of more than 2,300 K e pler c andidates is de- 
scribed in detail by iBatalha et al.l l|2012r) . The appendix 
table that contains the characteristics of the Kepler can- 
didates was extracted from the NASA Exoplanet Archivqj. 
We perform a similar calculation for At as described in the 
previous section. However, this time we calculate the differ- 
ence between fcirc and the duration measurement provided 
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Figure 1. The eccentricity distribution of the known radial velocity planets (top- left) and the calculated transit duration difference 
(circular vs eccentric) for b = 0.0 (top-right), 0.5 (bottom- left), and 0.8 (bottom-right). See ^for details. 
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Figure 2. The calculated transit duration difference (circular vs measured) for the Kepler candidates (left). The cumulative histograms 
for the Kepler candidates and the three radial velocity planets shown in Figurc[T]show a close match between the distributions, quantified 
by the K-S test described in Section [4] 
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by the candidates table, ikcpior- We do not use the provided 
h values since they are based on a circular orbit assumption 
from the measured transit duration and the stellar radii. 
We thus make no assumption on the value of b. We require 
Ri, > 0.7Rq to remove candidates for which the uncertainty 
in the stellar radius determination contributes significantly 
to the uncertainty in the transit duration. We also only in- 
clude candidates for which Rp > 8i?® to limit the sample to 
giant planets, similar to the RV sample. We show At versus 
a for the 176 resulting candidates in the left panel of Figure 

m 

To compare this distribution to its equivalents in Sec- 
tion [3] we perform a null hypothesis Kolmogorov-Smirnov 
(K-S) test to assess the statistical significance of their sim- 
ilarities. We binned the data by At into 40 equal bins of 
0.5 to collapse the data into 1-dimensional samples. The 
right panel of Figure [2] shows the normalized cumulative 
histograms for each of the four distributions: the three K-S 
tests compare the Kepler candidates to the RV planets with 
assumed b = 0.0 (Test 1), b = 0.5 (Test 2), and b = 0.8 
(Test 3). Test 1 produces a K-S statistic of _D = 0.05 that 
indicates a 100% probability that these data are consistent 
with being drawn from the same distribution (the null hy- 
pothesis). Test 2 produces a similar result of D — 0.075, 
also a probability of 100%. Test 3 results in D = 0.2 which 
is equivalent to a probability of 36%. This result can readily 
be seen in the right panel of Figure [2] where the 6 — 0.0 and 
6 — 0.5 cases are almost indistinguishable from the Kepler 
candidates, but the 6 — 0.8 case is clearly discrepant. As 
mentioned in 331 the small range of values for i for transits 
results in a uniform distribution of impact parameters with 
a mean of & = 0.5. The statistical congruence in the K-S 
test implies that the Kepler Mission is indeed recovering the 
eccentricity distribution of the RV planets. 

A criticism that may be levelled at this methodology is 
that the outcome of the statistical test depends upon the 
manner in which the data is binned. To determine the ro- 
bustness of our results, we used both half and double the 
number of bins to change the resolution of the sampling. 
For half the number of bins, we obtain D — 0.1 (100%), 
D = 0.1 (100%), and D = 0.2 (77%) for Tests 1, 2, and 3 
respectively as described above. If we then double the num- 
ber of bins, the results are D = 0.05 (100%), D = 0.0375 
(100%), and D = 0.1625 (22%) for Tests 1, 2, and 3 respec- 
tively. Clearly the results for Tests 1 and 2 are consistent 
with previous results and the results for Test 3 retain their 
discrepancies though with a variety of values. As described 
earlier. Test 3 (6 = 0.8) is the least relevant of the results 
since the mean impact parameter is b = 0.5. 

Finally, we investigate a sample of the outliers with 
particularly large deviations from the circular model (At > 
10 hours). These candidates are shown in Table [T] Since the 
Kepler data frequently do not contain any secondary eclipse, 
e and co are unknown. We calculate transit duration tccc as a 
function of e and lj via Equation [l] We then produce a grid 
of jtcirc — ikcpicrl/licirc — iccc | for all valucs of 6 and oj. Loca- 
tions where the grid values are approximately equal to 1 are 
possible solutions for which the measured transit duration 
in the Kepler candidate catalog is consistent with certain 
values of e and ui. 

An example of this is shown in Figure [3] where we 
present results of the above calculations as an intensity map 



Table 1. Minimum eccentricities for selected candidates. 



KOI 



Period 

(days) 



''kcplcr 

(hours) 



At 
(hours) 



44.01 


66.47 


19.74 


12.2 


0.74 


211.01 


372.11 


4.81 


10.5 


0.82 


625.01 


38.14 


4.24 


10.7 


0.85 


682.01 


562.14 


9.49 


10.8 


0.64 


1230.01 


165.72 


27.26 


11.8 


0.34 


1894.01 


5.29 


8.80 


14.4 


0.75 


2133.01 


6.25 


11.26 


14.1 


0.67 


2481.01 


33.85 


14.95 


16.8 


0.64 




Argument of Periastron (°) 

Figure 3. An intensity map for KOI 1230.01 that shows the result 
of dividing At using tkopicr by At using tccc- Thus, a value of 1 
(peak intensity) corresponds to the best solution (Q. 



for KOI 1230.01. In order to be compatible with the Kepler 
measured duration, the eccentricity of the planet must be at 
least 0.34. This process is repeated for each of the candidates 
in Table [l] in which we report the minimum required eccen- 
tricity emin for each candidate. It is worth noting, however, 
that these minimum eccentricities are not singular values 
but rather distributions, as can be seen by the gray-scale in 
Figure |3l The uncertainties depend highly upon the various 
random errors in the measured values of the Kepler candi- 
dates catalogue, including i. For example, the stellar radius 
of KOI 2481.01 would need to be ~ 45% of the catalogue 
value in order for it to be in a circular orbit and the duration 
discrepancy to be reduced to zero. 

Further of interest in Table [l] are the relatively short- 
period planets KOI 1894.01 and KOI 2133.01. One normally 
expects a transit duration of several hours for period such 
as these. However, the values of ikepier and At shown in this 
table imply a tcirc larger than 20 hours! This does not ap- 
pear to make sense until one considers the stellar radius. 
Note from Equation [l] that, for an edge-on orbit and small 
Rp, the transit duration scales linearly with the size of the 
star. For these two candidates, the stellar radii are 8.6 and 
9.3 solar radii respectively thus resulting in a large At and 
a significant eccentricity required to be consistent with ob- 
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2 4 6 

Minimum Planet Radius (Earth radii) 

Figure 4. The mean At for the Kepler candidates as a function 
of minimum planetary radius included in the sample (JJSj. 



servations. Note, however, that we have assumed b = 0. As 
one increases the impact parameter, the predicted transit 
duration will decrease and thus become closer to its mea- 
sured value. Results for individual cases extracted from the 
global distribution, such as those in Table [T] must therefore 
be treated with caution. 



5 PLANET SIZE CORRELATION 

We mentioned in Section [4] that the analysis of the Kepler 
objects included only candidates for which Rp > 8-R®. Here 
we perform a separate study by repeating the calculations 
of At for the Kepler candidates for a range of planetary 
radii. We allow the candidate sample to include all radii 
larger than lR(g to 8-R9 and calculate the mean of the At 
distribution in each case. We show our results in Figure [l] 

An interpretation of this figure is that the eccentric- 
ity distribution of exoplanets remains relatively flat until 
we probe below planets the size of Neptune. At that point 
the eccentricity distribution of the orbits becomes rapidly 
and significantly more circular. This is not unexpected since 
we understand from the solar system that the inefficiency 
of tidal dissipation (the quality factor Q) is much larger 
for h igh-mass than for low-mass planets IjGoldreich fc Soten 
119661 ). resulting in shorter tidal circularization time scales 
for smaller planets. One aspect of the Kepler candidate sam- 
ple that may influence this result is that they are dominated 
by planets at smaller semi-major axes since t hese have much 
larger t ransit probabilities. A s indi cated by iLissauer et al.l 
l|2011bl ') and ILissauer et al.l ((201^), multi-planet systems 
comprise a large proportion of the total Kepler candidate 
sample and these systems in particular are less prone to be 
false-positives. The findings that planet occur rence increases 
with decreasing planet mass (see for example [Howard et al.l 
I2OI0I ) then suggests that smaller planets find stable archi- 
tectures in systems with circular orbits and without large 
planets in eccentric orbits. This lends credence to two sce- 
narios: (1) core accretion forming terrestrial planets in cir- 
cular orbits, and (2) disk instability and capture scenario 
explaining the existence of giant planets in eccentric orbits. 

A potential alternative explanation for the dependence 



of At upon Rpianet IS a correlation between planetary radii 
and semi-major axes in the Kepler sample due to complete- 
ness. Smaller planets in closer orbits would be more likely 
to have had their orbits circularized leading to their domi- 
nation of the sample. We do not, however, see such an effect 
but instead find no correlation between planet size and pe- 
riod in the Kepler catalog for the selected range of orbital 
periods. Thus, we conclude that this effect is not the cause 
of the observed radius/eccentricity correlation. 



6 DISCUSSION 

There are various sources of potential systematic noise inher- 
ent in the data used to perform this analysis. For example, 
we have not taken the stellar limb-darkening into account 
when considering transit durations. However, since we only 
consider the total transit duration (first contact to last con- 
tact), this will be a negligible effect. 

For the Kepler candidates, the primary source of un- 
certainty arises from the stellar parameters that are used to 
derive many of the planetary candidate parameters. The pri- 
mary difficulties arise from the stellar radii whose precision 
is usually worse than ~10%, even when spectra are available. 
Our assumption is that these uncertainties are not signifi- 
cantly biased in one direction of the other and thus only add 
white noise to the overall statistical properties. One method 
to test this assumption is to consider multi-candidate sys- 
tems for which a change in the stellar radius will effect all 
calculated transit durations in a similar way. For example, 
KOI 157 (Kepler-11) has 6 detected planet candidates with 
measured transit durations that are shorter than predicted, 
and with an estimated host star radius of 1.06 Rq. Three 
of those discrepancies are quite small making them almost 
consistent with a circular orbit, as noted by Lissau er et al.l 
(|2011al ). Attempting to force circular orbits by reducing the 
stellar radius slightly makes 2 of the planets consistent with 
a circular orbit but leaves the other durations highly dis- 
crepant. A more global test of the assumption is the appli- 
cation of a range of uniform scaling factors to the stellar radii 
to try to produce transit durations consistent with circular 
orbits for the bulk of the distribution. These corrections did 
not change the distribution shown in Figure[2l leading to the 
conclusion that our calculated At values are not affected by 
systematically incorrect stellar radii. 

Considering the uncertainties in the radii of the Kepler 
host stars, our limits on the eccentricities of the specific Ke- 
pler candidates discussed in Section 2] should be treated in 
that context. Indeed, this is why we concentrate our com- 
ments on the global distribution of all the Kepler objects 
and their parameters, rather than on individual cases. Fur- 
thermore, the Kepler sample is only complete to ~ 0.5 AU 
with declining completeness beyond this to ~ 1.5 AU. Thus, 
the number of candidates beyond 0.5 AU is smaller but still 
sufficient for a valid comparison to be made. 

One aspect of the Kepler candidates that was not taken 
into account was the multiplicity of the systems. As sug- 
gested at the conclusion of the previous section, the multi- 
plicity may indeed play a significant role in stabilizing plan- 
ets in approximately circular orbits, particularly for those in 
the low mass/size regime. This is true of the radial velocity 
planets also, some of which are known to lie in multiple sys- 
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terns of super-Earth mass planets and with rela tively circu- 
lar or bits such as the system orbiting HD 10180 (jLovis et alj 

[mi. 

In Section [T] we men tioned the eccentricity bias found 
by IShen fc Turner! (|2008l ) due to low-amplitude signals in 
RV samples. This is a small but real effect which de- 
pends upon the sampling rate and has the consequence 
of under-estim a ting the number of near-circular orbits. 
IShen fc Turneil ([2003) develop a figure-of-merit and find 
that only ~ 10% of the planets considered in their sam- 
ple are affected by this bias. The samples studied here are 
two small to detect such an effect, but we mention it here as 
a consideration for future similar work for which the sample 
sizes and, more particularly, the period range explored have 
grown substantially. 

Finally, a minor impact on the stellar radii that should 
be noted is the one due to the relation betw e en pl anet 
frequency and stellar metallicity. Ijohnson et al.l (|2010r) ex- 
plored the mass-metallicity relationship for stars that har- 
bor planets and found a positive correlation of planet fre- 
quency with both stellar mass and metallicity, in accordance 
with the findings of iFisch er & Valenti (2005). This posi- 
tive correlation was also found empirically for M dwarfs by 
iTerrien et al.l l|2012f ). For a given stellar mass, a larger metal- 
licity leads to a smaller radius in order to reach hydrostatic 
equilibrium. The implication for this study is that many of 
the Kepler host stars will have relatively high metallicity 
leading to an over-estimated radius. However, this effect is 
at the level of a few percent and not expected to interfere 
with the results of this study. 



7 CONCLUSIONS 

By conducting a transit survey that is sensitive to long 
enough periods, it is expected that one will eventually repro- 
duce the eccentricity distribution found amongst radial ve- 
locity planets. This has not been possible until very recently, 
when the Kepler presented a large sample of long-period 
planets candidates, providing the incentiv e for this study 
and a similar one bv lPlavchan et al.l l|2012r ). For individual 
planets, the eccentricity may be d iscerned via a symmetry 
in the shape of ingress and egress l|Kipping||2008l ). This re- 
quires exquisite photometry and is highly sensitive to uj. We 
have shown here the consistency of the Kepler candidates' 
eccentricity distribution with their RV planets counterparts. 
The correlation of eccentricity with planet size is also an ex- 
pected result based upon the discoveries of small planets 
in multiple systems and indicates that there is an empiri- 
cal approach from which to both reverse-engineer formation 
scenarios and predict future stability patterns. 



planet Archive, which is operated by the California Institute 
of Technology, under contract with the National Aeronautics 
and Space Administration under the Exoplanet Exploration 
Program. 
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